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Abstract 

We study the effect of coherent propagation of two interacting particles in 
a disordered potential. The dependence of the enhancement factor for coher- 
ent localization length due to interaction is investigated numerically in the 
model of quantum chaos. The effect of interaction for two particles in many 

dimensions is also discussed. 

PACS numbers: 71.55.Jv, 72.10.Bg, 05.45.+b 



1 



Typeset using REVTgX 



I. INTRODUCTION 



The quantum localization of dynamical chaos has received a great deal of attention during 
the last years Q, 0. It has been understood that quantum interference effects lead to a 
suppression of diffusive spreading in the action space in spite of the chaotic dynamics of 
the correspondent classical model. An important consequence of this phenomenon is the 
exponential localization of quantum eigenfunctions over the unperturbed levels. A close 
correspondence has been established between this dynamical localization and the Anderson 
localization in a random potential for solid state systems || . One of the most studied models 
in this field is the kicked rotator model (KRM), which in the classical limit corresponds to the 
Chirikov standard map @ (CSM), a common paradigm of classical chaos. Although KRM 
seems to be at a first glance a pure mathematical model, it has however found important 
applications for real physical systems, as for example the process of microwave ionization of 
Rydberg atoms ||. Another useful property of KRM is that it can be studied very efficiently 
in numerical simulations allowing to investigate its properties in great details. 

Anyway KRM describes the one particle quantum dynamics and in many respects it is 
quite similar to one particle localization in a quasi one-dimensional random potential || . The 
last problem has been intensively studied in solid state physics and it is well understood from 
the theoretical point of view. On the contrary the case of interacting particles is much more 
complicated and a clear theoretical picture is still lacking (see for instance the recent Review 
0). Usually this problem is studied near the ground state and it is a common lore that a 
repulsive interaction would result in a even stronger localization f8|||. However, the recent 
investigation of two interacting particles in a quasi-one-dimensional random potential flO 



has shown that, even in the case of repulsive particles, interaction leads to an enhancement 
of localization length and a coherent propagation of two particles on a distance l c much 
larger than the one-particle localization length l\. According to [K| the enhancement factor 
is given by : 
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where U is the strength of on site interaction between two particles, V is the one-particle 
hopping element between nearest sites which determines the size of one-particle energy band 
and M is the number of transverse channels so that by itself l\ ~ M (M < In ( |1 . 1|) 
the intersite distance a is taken to be equal one and energy is taken near the center of the 
band so that fcp ~ 1/a = 1. Interaction U is considered to be less or comparable with V. 
The equation ( |1 . 1|) is valid in the regime when the enhancement factor is larger than one 
IJh > 1. 

The physical reasons according to which two particles propagates together, forming an 
effective pair, can be understood in the following way |T(J. Since the interaction couples 
only nearby sites (or on site) then the particles initially located on a distance r ^> l\ 
have only exponentially small effective coupling with each other due to the exponential 
decay of localized one-particle eigenstates. For such type of states two interacting particles 
(TIP) remain localized near their initial positions. These states are localized and form the 
majority of all states. However, there are other states in which TIP are initially close so 
that r <l\. As it follows from the previous case r » li, the initially close particles (r < li) 
cannot become separated on a distance much larger than l\ (otherwise we would contradict 
the previous case r ^> Therefore, they always remain on a distance r ~ li. On this 
distance interaction between particles is important. Qualitatively, the motion of one particle 
with respect to the other produces some kind of noise on it. This noise gives destruction 
of interference effects which have produced one-particle localization. The destruction of 
interference allows particles to propagate on a distance l c which is much larger than the 
distance between particles l\. This propagation can take place only for two particles since 
as soon as they become separated (r ^> li) they become localized. In some sense the way of 
TIP propagation in a random potential is similar to the method used by Miinchhausen to 
save himself from a swamp. 

The functional dependence in (|1 . 1| ) can be understood in the following way. The enhance- 



ment factor l c /h is proportional to the probability to mix one-particle states by interaction. 
In fact l c /l\ ~ Tp where T ~ U s 2 p is the rate of transition and p is the density of cou- 
pled states ||10|| . The coupling matrix element is U s ~ U/N 3 ^ 2 [[iTj] where N = Ml\ is 
the number of unperturbed components in one-particle eigenstate. The density of states 
is p = l/AE ~ N 2 /V so that the enhancement factor is l c /k ~ (UJAE) 2 ~ N(U/V) 2 
in agreement with ( |1 . 1|) . It is interesting to mention that such type of estimate is quite 
similar to the derivation of statistical enhancement for weak interaction and parity violation 



in neutron- nucleus reactions discussed by Sushkov and Flambaum |TT]. However in HIT 
even enhanced, the interaction remains relatively weak giving only small corrections while 
in our case the enhanced interaction U e ff ~ (Mli) 1 / 2 !/ is not small and it leads to significant 
changes in the properties of the system. 



Another way of derivation of (|1.1|) is based on the reduction of the TIP problem to 
a problem of superimposed band random matrices (SBRM) |H| . There, the interaction, 
even if repulsive, creates an effective thick wire along the diagonal ri\ = n 2 in the two- 
dimensional plane (ni,^) of indices corresponding to two particles. The effective width of 
the wire (effective number of transverse channels) is determined by the number of levels Ml\ 
coupled by interaction in one-particle basis. Outside of this width interaction is exponen- 
tially small and can be neglected at least in the first approximation. The large number of 
effective transverse channels M e ff = Ml\ leads to enhancement of localization length with 
the enhancement factor proportional to M e ff. 

All the above approaches were based on the assumption of statistical independence of 
transition matrix elements and eigenenergies in one-particle basis. This approximation seems 
to be reasonable due to randomness of potential and finite radius of interaction. However, it 
is very important to have a direct check and to verify the prediction ( |1.1| ). Some numerical 
checks were presented in [|TTJ. Here we present the results of a more detailed numerical 
investigation which we carried out for the model of interacting kicked rotators which had 



been discussed in flT0|]. We also present numerical results for a model with finite radius of 



interparticle interaction. 



The paper is constructed as follows. In section II we introduce the model and present the 
main results for on site interacting kicked rotators. The case of finite radius of interaction is 
discussed in section III. Conclusions and discussions of results are presented in section IV. 



II. THE "ON SITE" INTERACTION MODEL 



To investigate the effect of enhancement of localization length by interaction we used 



the model of two interacting kicked rotators introduced in | lQfl . The model represents two 
particles on a ring perturbed by kicks periodic in time. The evolution of the wave function 
ijj on one period of perturbation is described by the unitary operator (Floquet operator) : 

gt _ ^-iT{h\+h\)/2+iUS ni ,ri2 x ^-ik(cos 6»i+cos 2 ) ^2 X) 

with h\ i = —id I d6i 2- For U = we have two noninteracting kicked rotators which had been 
intensively studied during the last years 0, ||. The classical dynamics is chaotic and 
diffusive for the chaos parameter K = kT > 1 [Q. The diffusion rate is approximately D = 
n 2 /t = k 2 /2 for K ^> 1. Quantum interference effects lead to suppression of this diffusion 
for typical irrational values of T/Att and to exponential localization of eigenstates so that 
the averaged probability distribution over unperturbed levels decays as | ipn 1 2 ~ exp(— 2 | 
n — n | The localization length in the region of strong chaos is approximately given 

by l\ pd D w k 2 /2. Quasiclassical regime corresponds to k ^> 1,T <C l,kT = const and 
h > 1. 

For U 7^ interaction between particles is switched on. Using the Bessel expansion, eq. 
( |2.1D can be written as : 

g _ g— iT{h\+n^)/2+iUSni,n2 x j {fc)J ( ^)(— 2) ni+ra2_mi_m2 e* mi9l+ * m292 (2 2) 

mi ,m,2 

Since J n _ m (k) is exponentially small when \n — m\ > k, at each iteration of S many states 
(~ 2k) are coupled. Inter-particles interaction acts only when the two particles have the 
same momentum, namely when they occupy the same site on the momentum grid, on site 
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interaction, if we adopt the solid state terminology. Due to the presence in the exponent, 
the interaction U can only take values in the interval (0, 2ir). Due to interaction the two 
particles are able to propagate coherently on a distance l c much larger than the original one 
particle localization length 1%, as was anticipated in the introduction. Of course this can 
happen if they are initially started within a distance r < l x . Even if very close with the 



TIP problem in Id Anderson model ||10|| , our model has three different features. Indeed no 
randomness is here acting and the interaction is neither attractive, nor repulsive. In addition 
the perturbation couples many levels at each iteration (kick). 

The quantum dynamics was investigated in numerical simulations for symmetric config- 
urations with an effective number of unperturbed levels from 1000 to 2000. Antisymmetric 
configurations of two particles do not feel the on site interaction U and are not interesting. 
We iterated the quantum operator S starting from two particles initially at the same site 
for different parameters values. The spreading of the wave function in the 2-d space (ni, n 2 ) 
was studied through the second moments along the diagonal line n\ = n 2 : 

M*) = J((kl + M) 2 > t 

and across it 

*-(t) = <(M - \n 2 \f) t 

as a function of the iteration time t. In any investigated case a + was observed to saturate 
at an higher value than in absence on interaction, see Fig.l. On the other side <r_ keeps the 
same order of magnitude as l\ (as it should be in absence of interaction). This means that 
the localization length is strongly enhanced along the diagonal n\ = n 2 while it remains 
localized, with roughly the same localization length, across the diagonal. This is even more 
evident if one looks on the probability distribution P(rii,n 2 ) = \ip{ni, n 2 )\ 2 at a fixed time 
t ^> t*, where t* ~ l\ is the localization time, see for instance Fig.2. In this picture a local 
averaged distribution function is represented in the quarter of space ni, n 2 > 0, in a semilog 
plot. The channel of propagation along the diagonal n\ = n 2 is manifested in the contour 
lines drawn at the surface basis. 



From the distribution function important information can be extracted by computing 
the following distributions : 

P±(n ± )= £ m ni ,n 2 )\ 2 

|ni±ri2 \=n 

represented in Fig. 3 as a function of n± = \rii ± n 2 |/\/2 = n/^/2. These distributions give 
a measure of the "perturbed" localization lengths along (+) and across (— ) the principal 
diagonal. It is relatively easy to derive from them the respective localization lengths / ± 
by the usual best fitting procedure. Indeed the distributions are quite close to exponential 
curves P± ~ 2 exp(— 2n±// ± )// ± as can be inferred from Fig. 3. In the same picture we show 
the probability distribution in absence of interaction P° = 8n + exp(— 2 3 / 2 n + /7 1 )/7 1 2 with 
li = k 2 /2. This noninteracting distribution is quite similar to P_. The localization lengths 
are then plotted vs l\ = k 2 /2 in order to check the validity of Eq. ( |1 . 1|) . For sake of 
comparison the lines with power 1 and 2 are drawn. The dependence of coherent localization 
length l c = l + on one-particle length li can be satisfactory described by l + ~ 0.5/i 2 at U = 2 
while l~ m 1.5l\. However, the least square fit for the data of Fig.4 gives / ± ~ l\ a± with 
a + = 1.44 ± 0.29 and = 1.14 ± 0.07. We attribute the difference from the theoretical 
values 2 and 1 to the insufficiently large interval of variation of l\ (only 4 times). Further 
more detailed numerical investigations should be done to extract more accurate values of 
a±. Another interesting point following from the Figs. 1-4 is that for the same length li the 



coherent length l c is significantly larger than in TIP in Id Anderson model considered in |T0 
This can be seen by direct comparison of a + values. One of the reasons for this difference 
could be the different type of hopping in KRM where one kick couples many levels. 

To determine the numerical factor in the dependence of l c on both U and l\ one should 
also study the problem at small values of U 2. However, here for observation of the 
enhancement l c /h one should work at much larger values of l\ than we used in Figs.l- 
4. This requires a sharp increase of the basis and makes the numerical calculations too 
difficult. Therefore, to investigate the dependence on U we did the following. According 
to (|1 . 1|) we expect that it should exist a critical U cr given by C/crVn > C with C ~ 1. To 



check this we consider the same model but with random rotating phases, which means that 
T{n\ + n|)/2 in the first exponent is replaced by /(rti) + f{n 2 ) with /(n) = f(—n) being a 
random function in the interval (0,27r). 

In this way we can change configuration varying the random realization and obtaining 
results for the average behaviour. The results averaged over 10 realizations of disorder are 
presented in Fig. 5. The asymptotic value reached by the second moments a± = lim^oo a±(t) 
are plotted in units of the same value in absence of interaction (U = 0). Error bars are due 
to fluctuations in varying the random configuration. For small U, cr^° and crf° are both 
increasing up to double their value without interaction. When U > U cr full ( cr+) and open 
(er_) circles start to deviate one from each other thus indicating the presence of a sharp 
transition. In our case the transition starts at one particle localization length l x = 8 which 
approximately agrees with the observed critical value U cr ~ 0.3 and C ~ 1. We were not 
able to extract a more precise information on the dependence of l c on U due to the heaviness 
of numerical simulations. 

III. THE MODEL WITH FINITE RADIUS OF INTERACTION 

In this section we analyzed the effect of a finite range interaction on the dynamics. To 
be more precise we chose in (|2.1|) , instead of the former on site interaction U5 nitn2 a more 
general, finite radius interaction : 

Ur](n 1 ,n 2 )9(b - \n\ - n 2 \) 

where 6(x) is the usual step function which is zero for x < and one for x > 0. The phase 
7] is a random number in the interval (—1, 1) which depends only on n% if rii < n 2 and only 
on n 2 if n 2 < n\. It is quite clear that for b = it becomes the previous one with diagonal 
disorder. The diagonal disorder creates some difference from the model of section II since 
now the interaction depends not only from the difference n\ — n 2 . However, physically it is 
clear that diagonal disorder in interaction will not change too much the results. Indeed, the 
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main point is to have some coupling between two particles and the sign of interaction is not 
very important for the destruction of interference, since the one-particle random potential 
is already acting. Our numerical results confirm that disordered one site interaction gives 
qualitatively the same effects as for on site interaction U8 nuU2 . We usually investigated the 
cases with different interaction radius R = 2b + 1 and U = tt. 

Our main interest is to investigate the effect of interaction with finite radius R. From the 
theoretical point of view we can expect that for interaction radius R <l\ the equation ( |1 . 1|) 
is still valid since the particles are effectively coupled on a distance l\. However, for R > l\ 
the size of the effective thick wire on the lattice ni, ri2 is defined by R so we can expect that 
the enhancement factor will become larger l c /l\ ~ (R + li). This expression should remains 
valid up to values of R <C /12 where /12 is the 2d localization length for infinite radius R: 
ln/ 12 ~ l\. Indeed, for R Z 12 with the chosen type of interaction one should have the same 
localization length as in 2d. 

The results of our numerical simulations for finite interaction radius R are presented in 
Figs. 6-9. In Fig. 6 and Fig. 7 we show a + and cr_ for three different b values (b = 0,4, 16) 
which roughly agree with the above estimates. In agreement with the above picture the 
enhancement factor remains practically unchanged for R < l^. Only the case b = 16 has 
R > li and this produces a significant growth in a + (and even in <r_). The distributions 
P±(n), as defined in the previous section, are shown in Fig.8 and demonstrate a sharp 
increase of / + comparing to U = 0. In the same way we took the asymptotic values o~± 
reached by u±{t) at large time t and we plot in Fig.9 their square root as a function of the 
radius of interaction R (the values of have a similar behaviour). This figure confirms 
the above arguments that the enhancement starts to grow only for R > l\. However, it 
should be mentioned that the increase of R leads to a growth not only of o~ + but also of 
<7_. Indeed, a + and <r_ are growing in the same way: from the same Fig.9 one can see that 
the ratio a^/cr^ is approximately constant as a function of the interaction radius R. The 
physical explanation of this similar growth is quite simple: the increase of R leads not only 
to the increase of coherent propagation length but also to an increase of the effective size of 
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the pair which becomes of the order of R ^> l\. Unfortunately, we were not able to study 
numerically the regime R ^> h (in our case the maximal ratio R/h ~ 2) and it was not 
possible to check the dependence l c ~ Rl\. 

IV. CONCLUSIONS AND DISCUSSIONS 

Above we presented the results of our numerical investigation about two interacting 
kicked rotators in the domain of quantum chaos. They clearly demonstrate that on site 
interaction between two rotators in momentum space leads to large enhancement for local- 
ization length comparing to noninteracting case (Figs. 1-3). The localization length for 
coherent propagation of two particles l c = l + is significantly larger than the distance be- 
tween them l~ « l\. The maximal ratio l + /l~ in our numerical simulations was near 10 
(Fig. 4) which justifies the fact of effective enhancement of localization length for coherent 
propagation of two particles. The direct check of the relation ( |1 . 1| ) shows that the coherent 
localization length l c = l + grows approximately as l + ~ h 2 but more detailed numerical 
calculations are necessary to have a more accurate check of the power (see also discussion 
below) . 

Another part of our investigations was devoted to the effects of a final radius of interaction 
R between particles. They definitely show that for R < h the enhancement is not sensitive 
to the value of R (Fig. 9). The physical reason is quite clear. Indeed on site interaction 
couples one-particle states in a radius of h and therefore interaction with R < h does not 
give significant changes. For R 3> h the enhancement factor starts to grow with R. One 
can expect that in the regime i? 3> h the radius R will play the role of number of coupled 
states Ml\ = -R in an effective thick wire so that l c /h ~ R- Of course, this growth can 
continue only up to R < l±2 where lu is one-particle localization length in two dimensions 
and ln/12 ~ h ^> 1. While our results definitely show the increase of enhancement with 
R the power of growth is around 0.25 and is significantly less than 1. We attribute this 
difference to the fact that the ratio R/h was not big enough (in Fig. 9 R/h < 2.1) and the 



10 



asymptotic regime was not yet reached. The further increase of R is quite difficult since l c 
becomes comparable with the size of the basis. 

In general our results confirm the relation but a more detailed verification of this 
equation is still desirable. 

Let us now discuss in more details the different consequences of the result First 
we start from different dimensions d. For d = 2 the length l\ in ( [1 . 1|) should be understood 
as one-particle localization length in 2 dimensions. The number of transverse channels M is 
approximately equal to l± so that finally l c ~ /i 3 . For dimension d = 3 an interesting situation 
appears below Anderson transition for one particle [0. Indeed, it is possible to realize a 



random potential in which all one-particle eigenstates are localized for the hopping strength 
V < V c (a shift of mobility edge by interaction is not a very interesting case). As a typical 
example let us consider the Lloyd model with diagonal disorder i? ni n2jn3 = tan0 ni n2jn3 and 
hopping l^ona cubic lattice, where 4> ni ,n 2 ,n 3 are random phases homogeneously distributed 
in the interval [0, 7r]. In this case V c ~ 0.2 and below this value all states are localized. For 
two interacting particles in such random potential the effective strength of hopping for a pair 
will be strongly enhanced V e ff ~ y/NU. Here U is on site (or nearby site) interaction and 
iV ~ /i 3 is the effective number of states coupled by interaction. Since l\ can be quite large 
near (but below) one-particle transition point V c then two particles, even if characterized 
by repulsive interaction, can be delocalized when all one-particle states are exponentially 
localized |12| . Another way to see this effect is to say that the pair feels the disorder averaged 
over the size of the pair li which gives a strong effective decrease of disorder. Since in 3d 
derealization takes place for V e ff > V c generally there is no requirement to have l\ ^> 1 and 
it is not necessary to take V very close to V c . The condition V e ff > V c gives the boundary 
of pair derealization \Jl\^ 2 jV > 1. 

The appearance of derealization for a pair in 3d leads to quite interesting properties of 
energy spectrum. Indeed, for particles located on a distance r 3> li from each other the 
effective interaction is exponentially small ( ~ exp(— 2r/l\)) due to the small overlapping of 
one-particle states. Therefore, such states remain localized while the derealization will take 
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place only for the states with interparticle distance r < l\. Since the localized states with 
r ^> l\ form an everywhere dense spectrum this would mean that the continuous spectrum, 
corresponding to a delocalized pair, is embedded into the pure point spectrum of almost 
noninteracting one-particle states. 

Generally speaking such kind of spectrum is unstable with respect to small coupling 
between quasi-degenerate levels. In the present case the coupling is exponentially small 
but nevertheless it can change in principle the structure of the spectrum. The physical 
reason of such possible change can be understood in the following way. The delocalized pair 
propagates in a random potential which acts as some effective noise. This can increase the 
size of the pair even if the matrix elements for transitions with r = n_ ^> i x are exponentially 
small. Due to this noise the size of the pair will grow in time. The rate of growth can be 
estimated as D_ = n_ 2 /t ~ l\ exp(— 2n_/li). This gives a logarithmically slow growth of 
the pair size n_ ~ (^i/2) Int. At the moment it is not quite clear what will be the effect 
of the pair size growth on pair propagation in n + . At minimum, the displacement of the 
pair should become slower than diffusive n + 2 « {n\ — n 2 ) 2 ~ t/\nt. However, it is quite 
possible that sticking at n_ 3> l\ will produce a more significant effect on the growth of n + 
since in the region n_ ^> li the matrix elements for transitions in n + are also exponentially 
small. It is interesting to note that even in the case of strong attraction between particles 
the coupled state should be destroyed during the propagation in a random potential. Indeed, 
during the displacement of the pair disorder leads to transitions from the coupled state to 
continuum leading to the destruction of pair. Usually, the destruction rate is proportional to 
the squared amplitude of disorder and this can make the life time of coupled state relatively 
short. Contrary to this case the effective life time of a pair of repulsive particles discussed 
above can be much larger since n_ grows only logarithmically with time. In some sense 
the interference creates exponentially high barriers which effectively push particles to stay 
together. In quasi-one-dimensional case with l\ ^> 1 the effects of slow pair size growth can 
also lead to the appearance of logarithmic corrections in the expression of the enhancement 



factor in (|1 . 1|) . For example, we expect l c /h ~ /i/(ln/i) l/ with v ~ 1 for M ~ 1. The 
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effects of TIP in 3d systems below Anderson transition when all one-particle eigenstates are 
localized are quite interesting and at present, we try to study them in numerical simulations 
with effective 3d models ||13|| . Recently, an interesting approach to the TIP problem in 
(i-dimensions was introduced in |I4| . 

Up to now we discussed the effects of interaction only for two particles. However, for 
solid state systems the natural question is what will happens for a finite particles density 



p e . As it was discussed in |I| the above picture of TIP can be quite useful in the regime of 
small density !i C l/p c < l c . In this case the interaction is mainly reduced to interaction 
between two isolated particles. If all the particles are separated from each other by a distance 
L ~ l/p e ^> li then the interaction is exponentially small, all particles are localized and the 
current through such sample is exponentially small. However, it is possible to have another 
type of configuration when the particles are distributed by pairs of size 1%. In this case pairs 
can easily propagate on a distance l c 3> L ~ l/p e ^> l\. Collisions of pairs will go in a 
random way and will destroy interference effects for a pair. These collisions will lead to 
derealization and appearance of finite conductivity in an infinite system. It is interesting 
to note that it is enough to have only one pair when all other particles are well separated 
by the distance L 3> h- Then the collisions will allow to transfer the charge through 
the whole sample. However the above consideration, based on (|1 . 1|) and being correct for 
particle energy at the center of the band (E ~ V), should be applied more accurately for low 



energies near the ground state. Indeed, as it was discussed in [|1(| at low energies one should 



consider a transition from a lattice to a continuous system in which the enhancement factor 
should be proportional to l c jl\ ~ [kpl\)M since kpl\ determines the number of independent 
components in a localized state (for M 3> 1 the factor M should be replaced by kpcit where 
a t is the transverse width of the sample). Near the ground state p ~ kp and it seems that 
condition li <C l/p e ~ 1/kp implies that the enhancement does not work at low energies. 
Due to that at small densities there is no formal contradiction with the results || according 
to which repulsive interaction reduces the localization length near the ground state. To 
have a better understanding of the situation at small p e a more exact analysis should be 
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carried out to obtain a more precise expression for l c in the continuous limit. In principle, 
the average difference of energies for two repulsing particles (kp ~ 1/a — 1) on a distance 
r ^> h (Eqo) and r < l\ (-E/J is of the order of \Eoo — E^\ ~ U jl\ and is not very large 
for large l\. In fact this difference is even less than the amplitude of disorder W (we take 
the case of Id Anderson model discussed in [TTJ with diagonal disorder in the interval ±W 



where near the center of the band l\ ~ 25(V/W) 2 ). It is possible that for investigation of 
continuous limit <C 1 at low energy the approach used in |15|| for two particles with 
strong attraction can be useful after some extension. 

The most interesting case with density p e ~ 1 formally cannot be analysed on the 
basis of the result (|1.1|) for TIP. However, it is possible to think that interaction between 
quasi-particles can be studied in the same way as for TIP and that at small density of 
quasi-particles p q with k F ~ 1/a = 1 (i x <C l/p q <C l c ~ ^i 2 ) conductivity will be not 
exponentially small for one-dimensional samples with a size l sam l c ^> l\- In 3d for a "gas" 
of quasi-particles the possible slow growth of the pair size should be less important since 
collisions between pairs give rise to destruction of interference and finite conductivity in the 
regime where all quasi-particles are localized. Due to existence of exact connection between 
localization in Id and Id disordered spin systems [Q it would be interesting to understand 
possible manifestations of the analog of two particles interaction for spin systems. 

Finally, let us briefly discuss the possibilities of application of the observed enhancement 
for explication of large persistent currents observed in the experiments with small metallic 
rings ||16|| . Formally the coherent localization length ( |1.1|) is strongly enhanced in presence 
of interaction. Nevertheless, the direct estimates for the model of interacting kicked rotators 
and numerical results (see Fig.l and |10j) clearly show that the diffusion rate on the time 



scale l\ <C t <C l c is not larger than the classical rate at t <C h- It follows that the time to 
cross a sample will be not decreased by interaction. However, the magnitude of persistent 
current depends not only on the diffusion rate but also on the density of levels which in 
principle can become very large for multi-particles systems. Therefore, the possibility of 
enhancement of persistent current due to interaction is still open and should be studied in 
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more details. 
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FIGURES 

FIG. 1. Dependence of second moments on time in model (1.1) with k = 7, K = kT = 5, 
U = 2; upper curve is <r+, lower is er_. At t = both particles are at m = ri2 = 0, basis is 
-800 < n < 800. For U = <r+(t) » 600 for large t. 

FIG. 2. Probability distribution for two particles in the case of Fig.l at t = 8 x 10 4 . Different 
contours show different probability levels. 

FIG. 3. Probability distribution as a function of n± = 2~ 1 / 2 (ni ± 712) for the case of Fig.2: 
P + (n + ) (full line); P_(n_) (dashed); dotted line is the theoretical distribution P°(n+) for [7 = 0. 

FIG. 4. Dependence of localization length Z + (full circles) and l~ (open circles) on one-particle 
localization length l\ = k 2 /2 for K = 5, U = 2 and 4 < fe < 8. Full line shows dependence 
l c = l + oc l\ , dashed line marks l~ <xl\. 

FIG. 5. Dependence of enhancement for o± on U for the model (2.1) with random rotation 
phases (see section II), k = 4. Error bars are obtained from a± for 10 different realizations of 
disorder. 

FIG. 6. Dependence of a + on time for the model with finite interaction radius; k = 5.7, 
K = 5, U = it; R = 1 (full curve), 9 (dotted), 33 (dashed). Initial conditions are as in Fig.l, basis 
is —500 <n< 500. The noninteracting case U = has a + « 250 w l\ 2 (see Fig.3 in [10]). 

FIG. 7. Same as in Fig.6 but for cr_. 

FIG. 8. Probability distribution as a function of n± = 2~ 1 / 2 (ni ± ri2) for the case of Fig.6 
and R = 9: P + (ra+) (full line), P_(n_) (dashed line). P+(n+) (dotted line) is the theoretical 
distribution for U = 0. 

FIG. 9. Dependence of (full circles) and o°? (open circles) on interaction radius R; k = 5.7, 
K = 5, U = tt. 
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